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The existing Levy-Zhang approach to constructing the contribution to the resis- 
tivity of a metal of a magnetic domain wall is explored. The model equations are 
integrated analytically, giving a closed form expression for the resistivity when the 
current flows in the wall. The Boltzmann equation is solved analytically and the 
ratio of the spin up and spin down resistivities is calculated and its dependence on 
the strength of the Coulomb and exchange scattering potentials is elucidated. 
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I. INTRODUCTION 

Domain walls are examples of topological solitons in magnetism and they arise due to 
the competition between exchange and anisotropy energy. Domain wall motion by a spin 
polarized current has been gathering much interest recently mainly due to emerging device 
applications such as domain wall memory and domain wall logic devices. Winding number 
(vorticity), chirality and even skyrmion number are other degrees of freedom when consider- 
ing the magnetic domain wall and this is fascinating from the point of view of fundamentals 
as well as information storage considerations. Understanding the mechanisms by which a 
magnetic domain wall contributes to the resistivity of a metal, is a problem on equal footing 
with that of describing how a spin polarized current imparts torque to magnetization. When 
a conduction band electron fails to track the lattice magnetization when traversing a domain 
wall, an angle is subtended between the conduction band spin and the wall, which leads to 
a torque and, in the presence of impurity scattering, a measurable magnetoresistance. The 
relationship between domain wall motion spin transfer torque and domain wall magnetore- 
sistance was proposed by Tatara et al. 



2 



In this paper, we are interested in calculating analytically the explicit formula for the con- 
tribution of a domain wall to the resistivity in the diffusive limit, using the model equations 
of Levy and Zhang I2I. We wish to integrate this model, giving explicit formulae for the 
resistivity of a domain wall and use this formalism to calculate MR curves for systems in 
which domain walls nucleate in nanostructures by shape anisotropy. 



II. ADMIXTURE STATES AT A DOMAIN WALL 

We first begin with the simple picture of a 2 band ferromagnetic metal where the Fermi 
level lies in the Stoner split bands. We consider the Hamiltonian of a uniformly magnetized 
ferromagnet with the unit vector of magnetization aligned along the +z axis (cx.n = cr^), and 
so the starting SU(2) Hamiltonian takes the following form : 

Ho = ^^ + V{r) + Ja,, (1) 

where m* is the effective electron mass and V(T^) is the periodic crystal potential, taken to 
be invariant under SU(2) rotation and therefor this does not contribute to the spin scattering 
in the analysis which follows. We now write the Hamiltonian Hq in matrix form and look 
for eigenstates in the Hilbert space L2® Hs- 

where /? =h^/{2m*). We now transform Hq onto the basis {e*'^'^'''}, which assumes 
eigenvectors of the form $-^^=e*'^''' (ps^ where (ps is a two-component spinor. Writing 

hJ'''^' ° 1. (3) 

The eigenvalues of Equation [3] can be written as : \± =l3k'^ ± J where the ± signs refer to 
pure spin eigenstates. We now write the 2 component spinors for an unperturbed 2 band, 
exchange split ferromagnet. 



,(0) _ i_ I e^'^-^' \ (0) 1 / 



(4) 
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which describe pure spin states characterizing a two-band ferromagnet each with eigen en- 
ergies E| I = P k| ± J. We now turn our attention to the perturbation associated with 
a magnetic domain wall (DW), where the description of the conduction band spin goes 
beyond that of pure spins states. If the region of space over which the magnetization in 
a DW rotates is comparable the length scale of the Fermi wavelength l/kp, there will be 
an adiabatic 'mistracking' of the conduction spin with the lattice magnetization (see Fig. 
[1]). Nevertheless, as we will see, this scaling is treated as a perturbation in the Levy-Zhang 
approach. The perturbation parameter is proportional to VO.kp and this is assumed to 
be small for the perturbation expansion to converge. We write the Hamiltonian defined in 




FIG. 1: (Color Online) Schematic of a conduction band electron traversing a 180° Bloch domain 
wall and undergoing mistracking. 

equation [1] in the transformed basis of the domain wall Re$S°'' which gives : 

H = Rq ^HRq = Hq + Rg^ [H, Re] , (5) 

where Hq is now the unperturbed Hamiltonian of the magnetic system and Kg is the SU(2) 
rotation operator e~*i"'°". Indeed, in the basis of pure spin states, 6 is the polar angle 
of the magnetization unit vector n. We recognize that the perturbation potential can be 
written from the equation above as Vpert=Rg'^[H, Rg]. Now, Rg commutes with the J n.a 
and V(^)terms in Hq, so we are left with Vpert in the following form : 

Vp,rt = -PRg'[V\Rg] (6) 

= -/3Rg\\/''Re-RgV^) (7) 
= -P{R-^V'Rg-V'). (8) 



4 



Now, in order to evaluate the left hand side term of the above equation Rq ^V-Rg, we act on 
a trial wavefunction from the left as follows : 



{RfV.VRe)ij 
= Re'ViVRe^p + ReVi/j) 
= Rg\V^Re + {VRe)Vij + {VRe)Vij + ReV^i^) 
Rg\V^Reij + (Vi?e)V^ + {VRe)Vij + ReV^i^) 
= {Rg^V^Re + 2RfVReV + 



(9) 
(10) 

(11) 
(12) 
(13) 
(14) 



Inserting this into equation [HI we arrive at the following expression ; 



Vp,rt = -PiRe'^^Re - V=^)^ = -P{R^'\/^Re + 2R,\VRe)V + V')^ (15) 
We can can simply write VRg = — 'j^(n.(T)e~*^"'°" and V^Re = —^-^^{n.aYe~^^"'''^ - 



Recall that f3 = ff/2m* which is just a constant and that 6 is the angle of the magnetization. 
Recognizing that, for generators of SU(2) rotations, (n.a)^ = 1, which means that the first 
term in Vpert is diagonal and so does not mix spin states. Further, if the wall magnetiza- 
tion is assumed to be slowly varying in space with respect to the length scale defined by 
l/kp, we have V^e/{2VekF) < 1. This latter term may become important in DW profiles 
with vanishing V6 but finite (i.e. a stationary point in 6) which would occur in DW 
configurations with finite winding (n>l) or skyrmion number. As a first approximation, we 
retain the first order term in Vpert = -P'V6{n.a){—iV) and use the perturbation formalism 
outlined in Appendix A. We write the new eigenspinors in the rotated basis as : 



~iT:n.cr 



which, when substituted into Vpert now gives : 




(16) 




(17) 



(18) 
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For a Bloch wall, where the magnetization rotates in the yz plane, we now write the expansion 
coefficients for the first order corrections to the wavefunction (see Appendix A) , as follows : 

^ (fe(o)|-/3Vg.-^V(n.a)|nW) 

^(1) ^ ^ ry j^y nj± (21) 

(22) 

We find a similar expression for the Cj^'* mixing coefficient. It is important to note also that 
the unit vector along the magnetization can be written as n = (0, sin(^) sin(0), cos(^)) = 
(0, sech{x/X) sin(0), tanh(x/A)), for a Bloch wall in the +x direction with chirality 0. A is 
the equilibrium wall width (—^/lA/K) , A being the exchange stiffness and K is the magnetic 
anisotropy energy density). In this wall configuration, n.a has the components nya2 + 11^(73, 
where (Jj refers to the components of the Pauli spinors. Moreover, only the o"2 term yields a 
non-zero contribution to the mixing coefficient as it's elements are off diagonal. If the wall 
is set up to rotate in the xz plane, the coefficient C^^^ would be real. 



f 


^(0) _ ^(0) 

" k 




(tl- 


p(ye).{-iV9){h.a)\ 


i) 




p(0) p(0) 










-2 J + P{kl - A;f ) 



(1) _ /rf3r(i/3(e^(-'=T+fci)--V^.A;^ + ^e^(-*^T+fci)-'^)) 

~ /?(fc2-fc2)-2J ^^^^ 

(1) _ / ci'^(-i/5(e^(-'^i+^T)- v^.fc^ + ^e^^-^i+^t)- )) 

~ /5(fc2 - kf) + 2 J ^^^^ 

(25) 

In ferromagnetic metals, it is reasonable to assume that the kinetic energy splitting between 
the bands is much smaller than the exchange splitting , this condition is written as — 
1. It can further be assumed that ~ k^ , which is true to within an order of magnitude for 
most ferromagnets, and for the purposes of this calculation, the assumption is convenient 
in estabhshing the order of of magnitude of the effect. For spatially dependent V^, 
and n{pr).a this assumption would have to be relaxed and these terms will couple to the 
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scattering coefficients C^^^^ via a transformation on the basis kets {| t'), | J,')} via : 

(t' \zPVe.'k^h{-F').aV\ i), (26) 
a' I - i/3Ve.'k^h{r).aV\ T'), (27) 

which correspond to the Cj^'' and Cj^'* mixing coefficients, respectively. We now write 
the total wavefunction of the electron in terms of the adiabatically mixed two component 
spinors in the rotated basis, as follows : 



(28) 




(29) 



Using the approximations implemented by Levy and Zhang, we can write C^^^^ = and 
C|^'* = for a wall whose magnetization rotates along the x-axis in the adopted coordinate 
system. We now define the 'spin mistracking' parameter as : 



2J 



(30) 



Here is the local magnetization angle gradient, which can be taken to be locally constant. 
The normalization coefficients, can be written as follows ; N^j — = 1, we find that 

N(k^.)= 1 + (If)^- The value of can be taken to be locally constant over the lengthscale 
1/kp and for linearly varying magnetization profiles, the mistracking can be written as 



nx 
D 



(31) 

-D<x<D), originally 



AmJD 

which describes mistracking at the wall whose profile is - 9{x) 
considered by Levy and Zhang. We now write the corrections to the total wavefunction to 
first order as 



|*t( k , V)) = 77 



Rff 



Rf) 



^ fc^igi k f.'r' 
kp 

kp 

gi k |.T* 



(32) 
(33) 
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III. IMPURITY SCATTERING AT A DOMAIN WALL 



The domain wall itself does not necessarily give rise to inelastic scattering or to a mea- 
surable resistance. However, let us consider what happens when we consider an impurity 
potential to which the conduction spin is coupled via the coulomb and exchange interaction. 
The scattering potential is defined as follows : 

Vscatti^) = 5^ + JCT.M(ri)j - ^i) (34) 



This scattering potential has the following matrix elements in the basis as 
follows : 



>• (35) 
Using the basis defined by equation [331 we write down the matrix elements of the scattering 
potential : 



Rev I ^ I +ja.MiVi)Re 



5(T^ — 'r^i)d^~r^ 



(37) 



Note : I ^e^(k'-k).n|2 ^ ^ g.(k'-k).r.g^(k'-k).rj ^ ^ ^ ^ ^ g.(k'-k).(n-r,) ^ ^.^ ^^j^g^g 

i i,j i i^j 

Cj refers to the concentration of impurity scattering sites. This counting is an average over 
all of the impurity sites and each site is taken to be equivalent. Having established the 
matrix elements of the scattering potential, we write down the scattering rates based on 
Fermi's golden rule. After integration, the scattering matrix elements can be written in the 
following form : 

1 



|y-ni|2 



-1 2 



kk 



■N\K)N\k',) 



[v+j)k'^ - {v- ])k^ 



2^ 

k'p 



(38) 
(39) 



27r 



1 



kk' h 'N^{k,)m{k'^) 

271 



{v + aj) + 



n IP 



iv+j)k'^ - {v- i)k^ 



2^ 

k'p 



(40) 

(41) 
(42) 
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These scattering rates can be integrated over momentum space (k' coordinates) in order 
to find the total scattering hfetimes for momentum states within a spin channel and for 
momentum scattering which mixes the spin channels. This total scattering rate is defined 
as follows : 



(43) 



d^kN-^{k^)N~\k'^) 



{{v + j)K -{v- 3)KY^ + {{v + aj) + ^e{v - aj)) 

F F 



5(eka-ekv)- 



We expand the normalization constants N{kx) to second order ^ as follows ; N~'^{k.j.) ~ 
1 — {-^y + 0{{^Y) and we apply this approximation in order to evaluate the integral in 
Equation HH Recognizing that Equation HH has integrals of two types, we define these two 
types as follows ; 



/i = ^ / + Bf{l - + ..)5(^ic. - (44) 

/ d'k'^iC + D^ni - + ..)5(eic. - e,, J, (45) 

n J Kp Kp K,p 

(46) 

where we define the following constants within the integral. 

A = {v + j),B = -{v- j)k^, C={v + aj),D = k,{v - aj). 

The k space volume element is given in spherical polar coordinates as ci^k = 
k"^ sin(^^) cos{(j))dk and we write down Ii using this coordinate system, while expanding to 
second order in ^, as follows : 

1^ = 22^ f d'kk"^sm{e)d(f)dejj{Ak' sine' coscj)' + Bfx 
J kp 

(1 - 2klk'' sin^ e' cos^ <!>% + ..)5(ek. - 6,, J. 
/• t2 ^ (47) 

/ d' kk"^ sm{9)d(pdejj{A^k'^ sin^ 9' cos'^ + 2ABk' sine' cos(f)' + B^)x 
J kp 

(1 - 2klk'' sin2 e' cos2 + ..)5(ek. - e^' J- 



k ^^k 

We keep the approximation that the dimensionless mistracking is small such that -^rr^ <S 1, 

ftp 



and write the integral as ; 

27rc- /" / \ 

h = / d'kk'^5{ei^-ek,^>)d!0d'(f>'-^ (A^P sin^ cos^ 0' + sin^ 0' cos (/)' + sin 0'j 

(48) 

We now evaluate the integrals over k- space angle, as these are known analytically, as follows 



-'1 — 



J d'kk'^5{eka - ek'a')d'0d'(l)'^ (^A'k'^^ + 2AB7r.O + B\0^ (49) 



h 

which gives the result for Ii and this becomes, upon substitution for A ; 



k^^Hniv + j)2 

to order The integral of the term 5(ekcr — ekv) is constant with k-spacc polar angles 
as the band energy in the simplest case depends only on |kp and there is only a non-zero 
contribution to the integral for (5 (eta — e^'a') = 5{j3{k"^ — k'^) + (a — ex') J), in the case where 
we have k' =k and a' — a. We turn our attention now to the integral /2, which can be 
written as : 



h = '-f- j dh!h!H{e^„-e^,,,)d0'. 
d(f)' (C^ sin 0' + 2CL»E2!^Wfek! + k'^sinn'oos^d^'D^i^ ^ f^-^ _ 2(fc^^2 ^ 



(51) 



T 27rCi 

-'2 — 



J dk k ^5(ek(T — Gk'a')d(^ d(f) 
. (C^ sin 0' - sin^ O' cos^ cj^'^^^e + -g^''-^' ) (52) 

Evaluating these integrals over k space spherical polar angles, we now have : 



h = ^ J dk'k'^5{e^, - e^.„,)d0'd(t^ 
'c2.0 - c^f + '^CD=^ - - ik'^DC^^kl-^^ 

Integrating over k' and substituting in the definitions for C and D, we have : 



(53) 
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^ 3fc| J- ^"T^'^ + ^^^T^ ^^^^ 

We are now in a position to write down the total spin dependent scattering time as 

The equation above defines the momentum scattering time for the spin channel cr=± which 
refers to pure spins states. This is now used to solve the Blotzmann equation for the non- 
equilibrium distribution of electronic momentum which gives rise to the spin dependent 
diffusive current. 



IV. ANALYTICAL EXPRESSION FOR DW CONDUCTIVITY 

We start by finding the appropriate distribution function for the electrons in the metal, 
by writing down the first order solution to the Boltzmann equation as /*^k, t) = /'^(k — ^t). 
This is the distribution function for electrons in a field and the rate of change of this 
distribution function is given by : 

To first order in the electric field (E, with the convention e<0), we can write /°^(k) = /o(eko-), 
which simply corresponds to the Fermi-Dirac distribution and we now expand the last term 
in the Equation [56] above, as follows ; 

Vkr(ek) = :r^V6k. (57) 

= -hYlJiep-ei,^), (58) 
which gives us the field term in the distribution function rate equation : 

(%) =+eK-m^F-e^a)- (59) 

V / field 

We now turn our attention to the collision terms in the first order time derivative of the 
distribution function, and we begin with the spinless version : 



(%) = E Kk/(k')[l - /(k)] - l^kk'/(k)[l - /(k')]) . 

\ / coll , / 



(60) 
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where W^^'^ are the scattering rates (in units of energy per unit time). The first term in 
the equation above represents the 'scattering in' terms, while the second represents the 
'scattering out' terms and for elastic scattering, we have W^^' = W^^'k which arises due to 
the time-reversal symmetry which inelastic processes obey. 

(%) =E^kk'[/(k')-/(k)] (61) 

V / coll j^/ 

A J c^3kV,,4/(k') - /(k)] (62) 

Now we recast these equations in spin-dependent form, and now write the collision term, 
as follows ; 

(^) ,r ^ / K^^'f^''^^'^ " -^'^^^^^ + K^if'^i^) - n^)^) (63) 

and we invoke the steady-state condition : 



/ field \ / coll 

We now arrive at the appropriate Boltzmann equation for the spin dependent transport of 
electrons in the metal, in the most general sense : 

+ eK-ESie^ - ek.) = ^ / ^'k (l^^jr (k) - T (k')] + W^,'^[r{k) - /--(k')]) (65) 

+ ev^.E5(e^ - ek.) = J d'k\V^rf[r{k) - T (kO]5(ek. - e^v) + (66) 

K7\"iri^) - m^'M^^'T - 6k'_J (67) 

Now, for a two-band ferromagnet, we have the following relations for intra-spin-band and 
inter-spin-band scattering, respectively, as follows : 

eka - Ck'-a = ^ - ^''j + (69) 

Wc also expand the Boltzmann distribution about the Fermi energy and write the non- 
equilibrium distribution function for the electrons as follows : 

r(k) = /o(ek.) + n<^(k)5(e^ - ek.) (70) 



12 



47]-2 



J d'klV^rfln^imeF - - n-{^)S{eF - e^'J]6{e^^ - 6,, J + (71) 
|2[n-(k)5(6p - 6,v) - n-'^lkXe^ - 6,,„J]5(6k. - 6^,_J (72) 



lKTrK(k)-n-'^(k')]5(6^-6,,_J 



(73) 
(74) 



For the conductivity in the "current in wall" geometry, we solve the Boltzmann equation 
with the following distribution function : 

r (k) = r(k) - ev'^diep - ek.)r'^(k) (75) 

and write the total conductivity aciw = / — e/'^(k)vo-(i'^k of the system as 

<ycm = J d'k {-e\lE5{Ep - E^^)t\\,) - e\l5{Ep - E^{)t^{\,) + r(k)e(t;J + v^)) 

(76) 

The first two terms in the above equation are the spin dependant conductivities and 
a I respectively, while the last term is just the total conductivity of the two spin channels 
without spin relaxation. 

Takin the expressions for the total scattering rates from Equation HJl we set the scattering 
in momentum component to zero (the current is 'in wall') and write the the total spin 
relaxation time as : 



[r(k)] 



-1 _ 27rCi 



tip 



^sin 9' {v + ajY — {v + crj)'^k'^ sin^ 9' cos 



F 

2 



+[v + jyk"^ sin^ 9' cos^ d9' dcp' 



4,-27^ / dk'k"'5{ek^ - e;,.J ( Arriv + ajf - (v + cJjfk'^^l^ + (v + j^k'^^"" ^ 



3 k'jp 



l^e{A7r{v + ajf -{v + ajfPf^ + {v + j] 



3 k'j;, 



2 k^eHw 
3kp 



We now solve the Boltzmann equation, which is written as follows : 



-ev^E(5(e^-ek,) 



1 



in 



3 / 



r (k) - r (k') 



d^k'+ 



in 



3 / kk' 



r(k)-r'^(k' 



(77) 
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which, for the geometry under consideration, we write the solution to the Boltzmann 
equation as /"^(k) = /o(k) — eVyE5{eF — eko-)r'^(k), which is justified as integration over k' 
results in the vanishing of scattering in k components. We use this solution to write the 
conductivity of the wall with the current fiowing 'in wall' aciw- 

/ , " — — — — ^ (78) 

We now find that the total conductivity can be written as: 



W-^ / 42 42 

^ 4.42(.+,)2_(,+,)2l^ + (,+,-)2l^ 4<2(^_^.)2_(,_^.)2l^ + („+,)2l^ 



(79) 
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+ 



+ (^^+J 



^0) 



For the case for no n- vanishing k^^, where the current has a component perpendicular to 
the wall. Inserting the formulae for the spin scattering (relaxation) lifetimes from equation 
we get : 




6{Ef - ^kT)t^^k 



A;2 £-2 Stt 



(81) 



*2 



2 sin^ ^ sin^ <pd(pde5{EF - E],^)k'^dk 



(1- |sin2^cos2 



(82) 



sin^ 9 sin^ (Pd^d96{EF - Eki)k'^dk 
sm^9 sin^ ^d(pd96{EF - E^^^)k'^dk 



[iv+j] 



[V 



jy sin^ 9 cos^ 



sin^ ^ sin^ (j)d(f)d96{EF - Eki)k^dk 



[iv+j] 



jy sin^ 9 cos^ 



(83) 



As as result of the integration over k, we find that the expressions for the spin dependent 
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conductivity in the presence of a domain wall are as follows : 

= ^4^^ y (i; + j7(l-2sin2^cos2 0) ^ ^ 

—e'^h2'7iCi 2 3 /" sin^ ^ sin^ . 

^-t = J (^+^-)2l_4sin2^cOs2 0(^^-j)' ^ 

(86) 

Thus, the contribution to the resistivity is positive, as we have aciw = ctq + a-^ + ai and 
p = cr^^. The integral over angles in k space can be evaluated above to give the result: 

C^T = ^4^:^77 77711.6071, (87) 

while the cannot be integrated directly without prior knowledge of v and j. The current 
perpendicular to wall geometry (k±V^) can be solved. Let us now turn to the problem of 
the ratio of the conductivity of the spin channels, which is given by a =a^/ai, which we can 
write using Equations EH] as : 

sin^ sin^ 



J (t; + j7(l-2sin^cosV; , . 

sm t^sm 



(v + iY — 4 sin^ cos^ (piv — jY 
These integrals can be evaluated analytically, as follows : 

f-coseM - l-tan~\- ^^^ )y/cos{2e)cosec^{e) sin' e)\\l^ 
_2 2 V cos^ 6 



(pcosech^O tanh ^(j^ + 2vi + v"^) tan0 (gg) 

4(j - Vf - 4j3^2 _ 4j^3 + _ 2j4 COS 26 + 4j2t;2 coS 26 - 2V^ COS 26 

Sm U , ^= n 

40'--")^ ^i''-4j3D-10i2i,2_4jt,3+t,4_2j4 COS 20+4j2t,2 COS 26»+4j2D2 cos 26»-2i;4 cos 26* ^ 

Evaluating the right hand side of the above equation, we arrive at the simple relation : 

/•TT 

/ iv - ]f-nsixi6d6 , ..2 



/ 7rsin^d^(T;+j)' ^'''^•^ 



This establishes the dependency of the spin dependent conductivity asymmetry parameter 
a on the strength of the impurity potential v and the exchange coupling (j) to the im- 
purity. When there is an asymmetry in the impurity and exchange strengths (v^j), we 
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asymptotically approach the completely unpolarized current — ^0. Similarly, when the 
exchange potential strength vanishes (j — >0), we also tends towards the unpolarized case. 
On the contrary, when j=v, we now have a completely spin polarized case. The results of 
this calculation are plotted in Figure [2] 

CJ 

1 , 




FIG. 2: (Color Online) Plot of a Vs. the ratio of the impurity coulomb scattering potential and 
the exchange coupling to the impurity. 
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APPENDIX A: PERTURBATION EXPANSION TO SECOND ORDER 

We begin with a Hamiltonian which characterizes the ground state wavefunctions of the 
the total wavefunction ^) of the system . The Hamiltonian eigenvalue problem is written 
as follows : 

H\'^) = E\^) (Al) 

where |\E') is now written as an expansion of the ground states kets, which is linear super- 
position of orthogonal functions |n*^°^), lA;*-^-*), |m*^^)),... as follows 

(Ho+V){\n'')X\k^) + X^\m^) + ...) = {E^''^ + XE^^^ +E^^h^ + ....) {\n^)X\k^) + X^\m^) + ...) (A2) 

where the superindex refers to the order of the expansion. Since they are coefficients of 
an order-n polynomial, they are linearly independent. 



(0) 



n 



(0) 



FIG. 3: (Color Online) Spectrum of energy levels in a degenerate system. 
Taking the 0(1) equation from IA2l above, we have : 
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0(1) : Ho\n') = £;o|n°) (A3) 
^ (A4) 
{n^\Ho\n^) = Eo{n^\n^) (A5) 
0(A) : XHo\k^) + XV\n^) = XEo\k^) + XE^\n°) (A6) 
0(A2) : A'//o|rn') + AV|A;^) = A'£;i|A;^) + A'£;>°) (A7) 

(A8) 



We now arrive at the expansion coefficients for the first order correction to the wavefunctions 
, by multiplying the 0(A) equation above by the ket and we arrive at the following 
equation : 

A(n°|i/o|^') + A(n°|\/|n°) = XE^^\n^''^\k^^^) + XE^^\n^''^\n^°^) (A9) 
V n 7^ k, we have : 

(n(°)|l^|n(°)) = eW, (AlO) 

which is the first-order correction to the total energy. Taking the 0(A) equation, we 
multiply across by the ket \k^^^) armed with the decomposition of \k^^^) onto the vector 
space Ik^^'^), which we can write as |A;*^^)) = Ofc„ \k^'^'^). This can be written more 
succinctly in the outer product notation = \k^^^){n^^^\n^^^). We now write the 0(A) 

equation as : 

A(|A;(°)|i/o||A;(')) + A(|A;(o)|V|n(°)) = AEW(A;(o)|A;«) + AE«(A;(o)|nW). (All) 

Now, V n 7^ k, wc seek the cocfiicicnts of the expansion of the first order wavefunction in 
ground state kets, as follows : 

5^(A;(°)|i/oOfc|A;(°)) + {k^''^\V\n^''^) = E^^\k^^^Ck\k^^^) + (A12) 

=> (A13) 
= (A14) 
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In order to find the second order perturbation expansion coefficients, we iterative this 
progress further by using the O (A^), which has the following form : 

i/o|m(2)) + = £;(°V^'^) + ^^'^|^^'^) + ^^^'V^"^) (A15) 

//oC^|m(°)) + VC,\k^'^) = £;(°)C„|m(°)) + + (A16) 

(El^^ - El^y)C^\m('^ = VC,\k^'^ + + i?i^)|n(°)) (A17) 

Next, we multiply across by the ket |m°) and make use of the fact that the each basis 
\k^),\mP) are all orthonormal sets, as follows : 

Cm{E^^'> - ) = Vrr^kCk + E^^^CkSmk + E^^^ S^n (Al8) 



Next we take the case, whereby nT^k^^m, which gives : 

m+n ~ mjtn nj^k l"^"* ~ ) l"^" ~ ^k ) 

In the above equation, we have substituted in the result for the first order expansion coeffi- 
cients Cfe. The second order corrections to the total wavefunction are now know in terms of 
the matrix elements of the perturbing potential V and the ground state eigenenergies. To 

(2) 

find the second order corrections to the energy En , take n=m7^k in the above equation, 
which gives : 

Ei^^ = -KnmCfe (A20) 




FIG. 4: (Color Online) Exchange split ferromagnetic ground state. 
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